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ABSTRACT 

As part of a European Pulsar Network (EPN) multi-telescope observing campaign, 
we performed simultaneous multi- frequency observations at 1.4, 4.9 and 8.4 GHz 
during July 2006 and quasi-simultaneous multi-frequency observations from Decem- 
ber 2006 until July 2007 at 2.7, 4.9, 8.4, 14.6 and 32 GHz, in order to obtain flux 
density measurements and spectral features of the 5.5-sec radio-emitting magnetar 
AXP J1810— 197. We monitored the spectral evolution of its pulse shape which con- 
sists of a main pulse (MP) and an interpulse (IP). We present the flux density spectrum 
of the average proflle and of the separate pulse components of this flrst-known radio- 
emitting transient anomalous X-ray pulsar. We observe a decrease of the flux density 
by a factor of 10 within 8 months and follow the disappearance of one of the two main 
components. Although the spectrum is generally flat, we observe large fluctuations of 
the spectral index with time. For that reason we have made some measurements of 
modulation indices for individual pulses in order to also investigate the origin of these 
fluctuations. 



1 INTRODUCTION 



Magnet ars, as first discussed by iDuncan Thompson' 
(fl992. V are considered to be slowly rotating neutron stars 
with spin periods of 5-12 seconds and a rapid spin-down. 
They exhibit extremely strong magnetic fields, typically 
> 10^^ G, the decay of which is believed to create the observ- 
able high X-ray and gamma-ray luminosities, often visible 
in bursts. 

The model of magnetars tends to fit two previously 
described distinct classes of objects, the soft gamma- 
ray repeaters (SGRs) and the Anomalous X-ray pul- 
sars (AXPs). Common properties of the members of 
both classes are their long ro tation periods (AXPs ; 
iMereghetti Stellal[T995l ) . (SGRs; iKouveliotou et al.lll999l) 
and the bursting nature of their emission (Kaspi et al. 
l2003l : lGavriiL Kaspi Woods 2002 ). They are usually radio 
quiet. However, the pulsed ra dio emission that h as been de- 
tected from AX P J1810-197 (ICamilo et al.ll2006l ) and AXP 
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2003. With a period of 5.54 seconds and a period derivative 
- 1.15 X 10"^^ ss-\ a magnetic field of - 2.6 x 10^^ Gauss 
is implied. Although the previous observations clearly iden- 
tified it as an AXP, the extreme variation in the X-ray flux 
also classified it as the first transient AXP. 



1E1547.0-5408 (ICamilo et al 



l2008l ) 



indicates that this is 



not always the case. 

The AX P XT E J1810-197 was discovered by 
Ilbrahim et al. in Rossi X-ray Timing Explorer 

(RXTE) data of the source SCR 1806-20, taken in January 



The detection of a radi o source coincident w ith the posi- 
tion of AXP J1810-197bv lHalpern et al.1 (|2005l ) first raised 
the possibility that this was the first radio emitting mag- 
netar. This possibility was confirmed with the detection 
of strong, narrow and highly variable radio pulses, with 
the same pul se period as determined at high-energies, by 
ICamilo et al.l ([2006). Further observations have shown that 
the emission is ^-^80-95 % polarised, mostly linear, but with 
a significant degree of circular polarisation at all observed 
frequencies (jKramer et al.ll2007l : ICamilo et al.ll2007al V One 
of the most remarkable features of the radio emission from 
this source is its fiat radio spectrum. Radio pulsars are 
normally characteri sed by the steepnes s of their spectrum 
(Maron et al. 2000: lLQhmer et al.l[2008l ). which makes their 
detection at frequencies above 30 GHz very difficult. How- 
ever, soon after the discovery of its radio emission, it be- 
came clear that AXP J1810— 197had a fiat radio spectrum 
S oc zy-°-^ (jCamilo et al.1 l2006l . l2007d ) and for a time be- 
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came the brightest neutron star at frequencies greater than 
about 40 GHz. Strong variations in the pulsed intensity and 
the profile phase we re visible in the initial observations 
(jCamilo et al.ll2007bh . although it also became evident that 
the average pulsed flux density was decreasing with time. 

In this paper we report on the results of simultaneous 
and quasi-simultaneous multi-frequency observations con- 
ducted at the radio frequencies of 1.41, 4.90 and 8.35 GHz 
during July 2006 and 2.64, 4.85, 8.35, 14.6 and 32 GHz from 
December 2006 to July 2007. We present the radio spectrum 
of the total radio flux density of AXP J1810— 197 and of its 
individual profile components and the results of power-law 
fits to this spectrum. We also consider the time variability 
of the fitted spectral index in light of the modulation of in- 
dividual pulses and the intra-day flux density variability. A 
full discussion of the individual pulse properties will follow 
in a separate paper. 



2 OBSERVATIONS 

The simultaneous observations were made using the 100- 
m radiotelescope of the Max-Planck Institute for Radioas- 
tronomy (MPIfR) at Effelsberg, Germany, the 76-m Lovell 
radiotelescope at Jodrell Bank observatory of the Univer- 
sity of Manchester, UK and the 94-m equivalent Wester- 
bork Synthesis Radio Telescope (WSRT) in the Netherlands. 
The quasi-simultaneous observations were made with the Ef- 
felsberg radiotelescope. In total there were 8 simultaneous 
multi-telescope multi-frequency sessions during July 2006 
and 10 quasi-simultaneous multi-frequency between Decem- 
ber 2006 and July 2007. For the latter sessions the new sub- 
reflector of the Effelsberg telescope was used. It was installed 
in October 2006, improving the sensitivity and resulting in 
flatter gain curves but also allowing fast receiver changes 
between secondary and primary focus. The integration time 
for every session depended on the observing frequency and 
the observational circumstances. In general, the time needed 
for 1.42, 2.64, 4.85 and 4.90 GHz was around 5-15 min, for 
8.35 GHz around 20 min and for 14.6 and 32 GHz around 25- 
40 min. Details of the observing sessions are summarised in 
Tables ffland [2 

2.1 Calibration Procedures 

The observing and calibration procedures for Effelsberg 
data were the same for the simultaneous and quasi- 
simultaneous session. We used the 2.64 GHz, 4.85 GHz, 
8.35 GHz, 14.6 GHz and 32 GHz coole d HEMT receivers 
that are installed in the secondary focus (|Karastergiou et aP 
l200lh . In order to calibrate the flux density of a pulsar at 
the Effelsberg radio telescope reliably, a noise diode signal is 
injected into the wavefront following the horn synchronously 
with the pulse period in the beginning of each measurement 
(Figure El). The power output of the diode is compared with 
the power that is received from the pulsar. The signal of 
the noise diode itself is frequently calibrated and monitored 
by observing known reference sources during regular point- 
ing observations. The sources that were used for calibration 
were 3C273, 3C274, 3C286 and NGC7027. They were ob- 
served at the beginning and the end of each session, in com- 
bination with checks on pointing and focus stability. The 



Table 1. Summary of simultaneous observing sessions in July 
2006. 



Date 


Session 


Telescope 


Frequency 
(GHz) 


BW 
(MHz) 


31/05/06 


1 


Lovell 


1.42 


32 






WSRT 


4.90 


160 






Effelsberg 


8.35 


1000 


10/07/06 


2 


Lovell 


1.42 


32 






Effelsberg 


8.35 


1000 


17/07/06 


3 


Lovell 


1.42 


32 






WSRT 


4.90 


160 






Effelsberg 


8.35 


1000 


21/07/06 


4 


Lovell 


1.42 


32 






WSRT 


4.90 


160 


22/07/06 


5 


Lovell 


1.42 


32 






WSRT 


4.90 


160 






Effelsberg 


8.35 


1000 






Effelsberg 


14.60 


2000 


23/07/06 


6 


Lovell 


1.42 


32 






WSRT 


4.90 


160 


26/07/06 


7 


Lovell 


1.42 


32 






Effelsberg 


4.85 


500 






WSRT 


4.90 


160 






Effelsberg 


8.35 


1000 






Effelsberg 


14.60 


2000 


28/07/06 


8 


Lovell 


1.42 


32 






Effelsberg 


4.85 


500 






WSRT 


4.90 


160 






Effelsberg 


8.35 


1000 






Effelsberg 


14.60 


2000 



whole procedure o f Eff elsberg calibration is described fur- 
ther in Angelakis (2007). We also monitored the quality of 
our observations by taking data for well known pulsars such 
as PSR B 1929+ 10 which is strong enough to be detected at 
all frequencies up to 43 GHz. The observed properties of this 
pulsar, extensively also studied at Effelsberg, were compared 
to archival data to confirm that the system was functioning 
correctly at all frequencies. The observing and calibration 
procedures for the other telescopes participating in our si- 
multaneous observations between May and July 2006 are de - 
scribed for each telescope in detail in .Kramer et al.l (|2007l ). 



3 DATA ANALYSIS & RESULTS 

Until August 2006, the integrated pulse profile of 
AXP J1810— 197 consisted of two major well-separated fea- 
tures ^which we discussed in some detail already in 
Kram er" et al.l (2007). Following the same convention, we re- 
fer to them as the main pulse (MP, right feature in Figure[2^) 
and interpulse (IP, left feature in Figure [2^). The MP is the 
more complex and wider of the two, with a varying width 
of about 0.7 seconds (or 45 deg in pulse longitude) while 
the IP is much narrower with a width of typically 0.2 sec- 
onds (or 12 deg longitude). In some cases, three or more dis- 
tinct sub-components were visible in the MP. The simpler 
IP was not always visible and was only strongly detected 
during parts of our observations, as indicated in Figure [1] 
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Table 2. Summary of quasi-simultaneous observing sessions from 
December 2006 to July 2007. 



Date 


Session 


Telescope 


Frequency 
(GHz) 


BW 
(MHz) 


09/12/06 


1 


Effelsberg 


2.64 


100 








4.85 


500 








o.oO 


1 nnn 

lUUU 








14.60 


2000 


ZD / iZ/ UO 


9 
Z 


rL/iielsberg 




1 nn 

lUU 








4.85 


500 








8.35 


1000 








14.60 


2000 


n/1 /n9 /n7 
u^/ uz / U ( 


q 


rL/iielsberg 




1 nn 

lUU 








4.85 


500 








8.35 


1000 








14.60 


2000 


/n9 /n7 
UO/Uz/U < 


/I 


Effelsberg 


/I QC^ 
4.00 


i^nn 








32.00 


2000 


12/02/07 


5 


Effelsberg 


4.85 


500 








14. ou 


9nnn 
zuuu 


17/02/07 


6 


Effelsberg 


4.85 


500 








8.35 


1000 








14. ou 


9nnn 
zuuu 








32.00 


2000 


18/02/07 


7 


Effeslberg 


4.85 


500 








8.35 


1000 








14.60 


2000 








32.00 


2000 


26/03/07 


8 


Effelsberg 


4.85 


500 








14.60 


2000 








32.00 


2000 


05/05/07 


9 


Effelsberg 


4.85 


500 








8.35 


1000 








32.00 


2000 


06/07/07 


10 


Effelsberg 


4.85 


500 








8.35 


1000 








14.60 


2000 



where we present the daily average flux densities of the two 
components as measured at 8.4 GHz as a function of time. 
The identification of the remaining visible pulse feature (see 
e.g. Figure [SId) with the MP is possible due to timing infor- 
mation obtained from regular monitoring observations with 
the Lovell telescope at the Jodrell Bank observatory (Lyne 
et al. in preparation). After the disappearance of the IP in 
summer 2006, both MP and IP were visible again simulta- 
neously only during one short session in May 2007. The IP 
remains undetected since. 

Overall, the flux density of the source has significantly 
decreased since its first detection at radio frequencies (see 
Figure [D , wh i ch is consistent with the earlier findings by 
ICamilo et al.l (|2007bl ) . Due to differences in the frequency 
range of the simultaneous and quasi-simultaneous observa- 
tions we discuss further results from these campaigns sepa- 
rately. 



■ Main Pulse 
A Interpulse 



54100 

MJD 



Figure 1. The average flux density for each observing session 
as measured at 8.35 GHz for the main pulse and the interpulse 
(when detected). 



3.1 Simultaneous Observations 

The flux densities measured during the simultaneous obser- 
vations and the spectral indices derived from power law fits, 
S (X ly^) to these data, are Usted in Tableland are shown in 
Figures [3] and IH For each observation we summed all indi- 
vidual pulses to obtain an integrated total power profile that 
was flux calibrated. The nominal conservative uncertainties 
for our flux density measurements are about 10% (see also 
Section 3.2). In addition to our flux densities, we also utilize 
published flux densities measu red at 1.42 and 4.8 GHz with 
the VLA (jCamilo et al.ll2007dl and at 1.4 and 4 8 GHz at the 
Mount Pleasant observatorv ("Hotan et alllioOTl ) (see Fig. [3)). 
Our flux densities observed for the MP at similar epochs for 
these frequencies are in good agreement with these values. 
This fact indicates that the effect of interstellar scintillation 
may be small due to the large observing band widths, but we 
will discuss the possible impact of the interstellar medium 
on our results in more detail later. Overall, a large decrease 
in the signal strength with time is observed. At the same 
time, the spectrum is generally flat with an average spectral 
index of a = -0.31 =b 0.06 for the MP and a = -0.41 ± 0.21 
for the IP (s e e Fig , p and ID, consistent with results by 
ICamilo et all (|2006l . [2007ch . However, around MJD 53938 
(2007 July 22) we observe a strong time variability, with the 
spectrum being initially steep and then flattening day by 
day for both MP and IP. 



3.2 Quasi-Simultaneous Observations 

All quasi-simultaneous observations were performed with Ef- 
felsberg alone. The short time needed to switch between the 
different receivers (i.e. ^ 30 — 60s between secondary focus 
receivers, and ^ 2 — 4min between secondary and primary 
receivers) made it possible to observe at a wide range of 
frequencies in a single observing session. The longest multi- 
wavelength observing session (without the calibration scans) 
of AXP J1810— 197 was about four hours, during which fre- 
quencies were cycled through repeatedly in order to detect 
any short-term variability. Using these data, we not only 
studied the long term variability of AXP J1810-197 , but 
also the medium-term intra-day flux density fluctuations as 
well as the modulation of individual pulses during each of 
these sessions. An example of such measurements is shown in 
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Table 3. Summary of flux densities and spectral indices for the 
main pulse. 



Date Frequency 

(GHz) (mJy) 



31/05/06 


1.42 


8.23 ±0.82 


±0.01 ±0.11 




4.90 


9.73 ±0.97 






8.35 


8.03 ±0.01 




10/07/06 


1.42 


6.90 ± 0.69 


-0.27 ± 0.06 




8.35 


4.28 ±0.03 




1 7 /n7 /nfi 


1/19 


D.UD zt U.Di 


— zlz U. 




A Qfl 


9 97 -1- n 9^^ 






O.OO 


9 _|_ n no 




21/07/06 


1.42 


2.14 ±0.21 


-0.33 ±0.11 




4.90 


1.42 ±0.14 




22/07/06 


1.42 


2.59 ±0.26 


-0.79 ±0.30 




8.35 


1.04 ±0.03 






14.60 


0.33 ±0.03 




23/07/06 


1.42 


3.24 ±0.32 


-0.55 ±0.11 




4.90 


1.61 ±0.16 




26/07/06 


1.42 


0.99 ±0.10 


-0.03 ±0.11 




4.85 


1.40 ±0.03 






4.90 


0.86 ±0.09 






8.35 


1.15 ±0.03 






14.60 


1.06 ±0.04 




28/07/06 


1.42 


1.32 ±0.13 


0.15 ±0.28 




4.85 


1.66 ±0.03 






4.90 


0.65 ±0.06 






8.35 


1.92 ±0.07 






14.60 


1.70 ±0.04 





Table 4. Summary of flux densities and spectral indices for the 
interpulse. 



Date Frequency 

(GHz) (mJy) 



10/07/06 


8.35 


0.86 ±0.02 




17/07/06 


1.42 


0.13 ±0.01 


±1.12 ±0.33 




4.90 


0.83 ±0.08 






8.35 


0.86 ±0.02 




21/07/06 


4.90 


0.44 ± 0.04 




22/07/06 


4.90 


0.58 ±0.06 


-1.61 ±0.08 




8.35 


0.27 ±0.02 






14.60 


0.10 ±0.04 




23/07/06 


4.90 


0.58 ±0.06 




26/07/06 


4.85 


0.81 ±0.02 


-1.29 ±0.26 




8.35 


0.51 ±0.02 






14.60 


0.19 ±0.03 




28/07/06 


1.42 


0.49 ±0.05 


0.15 ±0.28 




4.85 


0.42 ±0.02 






4.90 


0.28 ±0.03 






8.35 


0.85 ±0.04 






14.60 


0.58 ±0.03 





"5^0 1^60 1^^90 2T2O 2^50 slso 3^10 So 4^70 msec 

28/ 7 / 2006 
8.350 GHz 

,1 




.0^ 34.4 68.9 103.3 137.7 172.2 £06.6 241,1 £75.5 309.9 (Peg.) 

I 5?0 1^60 iSo 2T2O 2^50 slso 3^10 4S0 4^70 msec 



9 / 12/ 2006 
8.350 GHz 




"sS 1^60 1S0 2T2O 2^50 slso 3^10 So 4^70 msec 



5 / 5 / 2007 
8.350 GHz 




.0| 34.4 68.9 103.3 137.7 172.2 £06.6 £41.1 £75.5 309.9 (Peg.) 



Figure 2. Phase aligned integrated proflles of AXP J1810— 197in 
2006 July (top), 2006 December (middle) and 2007 May (bottom). 
The source becomes signiflcantly weaker while the IP eventually 
disappears. Only in 2007 May, both the MP and the IP are visible 
again. Note that each shown proflle contains the signal of the 
calibrating noise diode which was switched on synchronously with 
the pulse period for the flrst 50 phase bins. 



Fig [5] where we present the flux density as measured repeat- 
edly at 4.85 GHz over two consecutive days. The observed 
variation is consistent with a constant flux density over this 
period. The results of all our quasi-simultaneous observa- 
tions are summarized in Table [5] 

The average flux density values were then used to deter- 
mine the spectra and their indices. This was done for each 
session and the results are shown in Figure (6] The results 
are consistent with and extend those of the simultaneous 
measurements. We find again that the spectrum is generally 
flat with a mean spectral index of a = 0.00 ± 0.09 but it is 
also variable on a day-to-day basis with significant variations 
around its mean value. This is reflected also by the varia- 
tion of the flux densities measured at 8.35 GHz and above 
which appear to sometimes follow a common pattern that 
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Figure 3. (Top) Flux density spectrum of the MP as determined 
during our simultaneous observations. For plotting purposes the 
flux densities of each day were multiplied by a different factor 
to distinguish the dates more clearly. (Middle) Flux density and 
its variability as measured at various frequencies as a function 
of time. We also added data measured by Hotan et al. (2007) 
at 1.4 and 4.8 GHz (shown as crosses and stars, respectively) and 
Camilo et al. (2007c) at 1.42 and 4.8 GHz (shown as bars and exes, 
respectively). (Bottom) Derived spectral indices of the MP as a 
function of time. The average value is determined a = —0.31 ± 
0.06. In comparison, we show spectral indices derived by Camilo 
et al. (2007c) from VLA measurements as triangles. 



is anti-correlated with changes seen below that frequency. 
There is an indication of a general increase of a over the 
210 days covered by our observations, starting with initial 
values around a ~ —0.2 and reaching a ^ 1.2 near the end 
of the observing period. 
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4.90 GHz 
• 4.85 GHz 




1 


J. 




■ 1.42 GHz 










Figure 4. (Top) Flux density spectrum of the IP as determined 
during our simultaneous observations. For plotting purposes the 
flux densities of each day were multiplied by a different factor 
to distinguish the dates more clearly. (Middle) Flux density and 
its variability as measured at various frequencies as a function of 
time. (Bottom) Derived spectral indices of the IP as a function 
oft time. The average value is determined a = —0.41 ±0.21. 



3.3 Modulation Indices 

In order to characterise the observed flux density varia- 
tions and to obtain reliable estimates for the measurement 
uncertainties, we also studied the modulation of the sin- 
qle pulse flux densi ties (see e.g. IKramer et al.l (|2003l ) and 
IServlak et all (|2007l ll. We can use these results also to esti- 
mate the impact of the interstellar medium on our measure- 
ments, following similar studies such as that conducted for 
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54079.0! 



-20 20 40 60 80 100 120 140 160 180 200 

Time (min) 

Figure 5. Flux density measured at 4.85 GHz during an ob- 
serving session for two consecutive days (square and up triangle 
points). The first two measurements for each day correspond to 
a 5 minutes integration time while the third measurement corre- 
sponds to 15 minutes. 



Table 5. Summary of flux densities and spectral indices for the 
MP measured during our quasi-simultaneous observing sessions. 



Date 


Frequency 
(GHz) 


Sm 

(mJy) 


a 


09/12/06 


2.64 


0.39 ±0.05 


-0.35 ±0.11 




4.85 


0.40 ± 0.03 






8.35 


0.27 ±0.01 






14.60 


0.23 ±0.03 




26/12/06 


2.64 


0.56 ±0.09 


-0.67 ±0.17 




4.85 


0.25 ±0.05 






8.35 


0.24 ±0.05 






14.60 


0.16 ±0.01 




04/02/07 


2.64 


0.55 ±0.08 


-0.24 ±0.49 




4.85 


0.15 ±0.05 






8.35 


0.15 ±0.01 






14.60 


0.36 ±0.03 




06/02/07 


4.85 


0.20 ±0.05 


-0.18 ±0.25 




32.00 


0.14 ±0.06 




12/02/07 


4.85 


0.20 ±0.03 


0.11 ±0.15 




14.60 


0.22 ±0.02 




17/02/07 


4.85 


0.24 ±0.02 


±0.06 ±0.13 




8.35 


0.18 ±0.03 






14.60 


0.23 ±0.01 






32.00 


0.25 ±0.05 




18/02/07 


4.85 


0.24 ±0.03 


±0.03 ±0.08 




8.35 


0.20 ±0.03 






14.60 


0.21 ±0.02 






32.00 


0.25 ±0.06 




26/03/07 


4.85 


0.33 ±0.03 


-0.80 ±0.17 




14.60 


0.10 ±0.02 






32.00 


0.08 ± 0.04 




05/05/07 


4.85 


0.09 ±0.03 


±0.86 ±0.27 




14.60 


0.09 ±0.05 






32.00 


0.42 ±0.08 




06/07/07 


4.85 


0.15 ±0.04 


±1.20 ±0.44 




8.35 


0.19 ±0.03 






14.60 


0.55 ±0.02 






Frequency (GHz) 



-i — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — |— 
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540600 540900 541200 541500 541800 542100 542400 542700 543000 

MJD 

Figure 6. (Top) Flux density spectrum of the MP as determined 
during our quasi-simultaneous observations. For plotting purposes 
the flux densities for each day were multiplied by a different fac- 
tor to distinguish the dates more clearly. (Middle) Flux density 
and its variability as measured at various frequencies as a func- 
tion of time. (Bottom) Derived spectral indices of the MP as a 
function oft time. The nominal average value is determined as 
a = -0.003 ±0.089. 



pulsars bv lMalofeev et aT] (|l996h at very similar frequencies, 
i.e. 4.75 GHz and 10.55 GHz. 

We concentrate on the flux densities measured at 
4.85 GHz and 14.6 GHz during our quasi-simultaneous ses- 
sions. The flux densities at these frequencies were measured 
in every session, often at the beginning and the end of a 
session, providing a densely sampled data set that gives a 
good representation of the behaviour between widely spaced 
frequencies. 

As a comparison, we also studied data for our reference 
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Figure 7. Variation of the pulse-to-pulse modulation index of 
AXP J1810-197at 4.85 GHz (squares) and 14.6 GHz (triangle). 
At 14.6 GHz the indices were computed for time series covering 
25-30 min while the 4.85 GHz data cover 10-30 mins. 



source PSR B 1929+ 10 to rule out systematic effects, impact 
of weather, or instabilities in the receivers chain. We find, 
as expected, only minor variations in the flux density of 
this well k nown pulsar at high fr equencies, in line with the 
findings of iMalofeev et al. 

For each session, we calculated the pulse-to-pulse mod- 

- — Tsy^ — - < 
where S is the measured 



ulation index m according to 
S >^ (e.g. Kramer et al. 2003), 
flux density, (S) its mean value and as its standard devia- 
tion. We present the resul ts in FigurelT] As already expected 
from our discussion in Kr amer et al.l (j2007) and co nsisten t 
with e arly observations of AXP J1810— 197 bv lCamilo et al.l 
(|2006l ), the pulses are very highly modulated and variable 
for both frequencies. For the densely covered period of time 
around epoch MJD 54150 (04 February-26 March 2007) the 
data occasionally suggest a behaviour that is anti-correlated 
between 4.85 GHz and 14.6 GHz. Overall, the higher fre- 
quency shows greater modulation, but it is possible that 
weaker single pulses are more difficult to detect at these fre- 
quencies. Despite this large modulation of the single pulses, 
averaging over sufficient time as done during all our observa- 
tions essentially removes the variation for a given measure- 
ment value, leading to relative uncertainties in S of the order 
m/ ^Jn where n is the number of averaged pulses. This is con- 
sistent with our estimated flux density errors and conflrmed 
by the repeatability of our flux density measurements over 
consecutive days (see Fig [5}. We therefore expect the im- 
pact of slow- varying effects caused by refractive scintillation 
to be more important. We will discuss this in the following. 



4 DISCUSSION 

Our monitoring observations of AXP J1810— 197 and its flux 
density spectrum over the course of more than one year ap- 
pears to conflrm the conclusion derived already in all previ- 
ous studies of it, that AXP J 18 10— 197 behaves unlike most 
or even any known radio pulsar. Several observed properties 
support this fact. In our range of studied frequencies, the 
source exhibits a generally flat flux density spectrum, with 
values of spectral index co nsistent with those measured by 
ICamilo et all (|2006l . l20073 l . These flat spectra allowed us to 
observe the source up to very high frequencies, at least when 



compared to radio frequencies typically used for ordinary 
pulsar observations. An example is our detection of the mag- 
netar at 43 GHz in May 2007 with the Effelsberg radiotele- 
scope. This is only t he flfth neutron star detected at A7mm 
()Kramer et al.lll997h and the only one where single pulses 
could be observed. This is consistent with o bservations at 
the IRA M 30 m telescope at 88 and 144 GHz bv lCamilo et al.1 
(|2007d ). representing the highest radio frequency at which 
a neutron star has been detected. For normal pulsars, we 
often see breaks in their power law spectra at high fre- 
quen cies (jMaron et al.l I2OOOI : iKiiak. Gupta h Krzeszowskil 
I2OO7II or even occas sionally a turn-up at mm- wavelengths 
([Kramer e t al.|[l997l ). Here, in our range of frequencies, the 
spectrum of AXP J1810— 197 is well described by a single, 
flat power law, of the form Sv oc with a — 0.0 =b 0.5. 

In those cases where both the MP and IP were detected, 
their spectral indices were found to be differing, with the 
IP showing a signiflcantly greater variation. Even though 
the spectra appear to change in the same manner from day 
to day, becoming steeper or flatter together, the IP exhib- 
ited many more extreme positive and negative steep values. 
This is somewhat surprising since the MP spectral index 
reflects emission from a wider range of emission compo- 
nents, in comparison to an IP that is always observed as 
a rather simple emission feature. The variation of the spec- 
trum is roughly consistent with the model by Thompson 
(2008) where timescales of ^ 0.1 day may be explained by 
current- driven instabilities on the closed magnetic fleld lines. 
The observation that it affects the IP and MP in a differ- 
ent way is interesting and may support the idea of a rather 
twisted magnet osphere. Although the spectrum is generally 
flat, the spectral index shows signiflcant variation, with a 
slight trend of becoming more positive with time. 

While we have shown that the effect of the signifi- 
cant pulse-to-pulse modulation can be removed by aver- 
aging a sufficient number of pulses, another well known 
cause for flux density variations, and hence possible spectral 
index variations, is interstellar scintillations (ISS). Based 
on studies of t he IS S at 4.75 GHz and 10.55 GHz by 
iMalofeev et aD (119961), the relatively large dispersion mea- 
sure of DM = 178pccm~^ and the estimated distance 
of 3.3 kpc, suggests a critical scintillation frequency /c of 
~14GHz, well within our observing band at 14.6 GHz. At 
this frequency, we may therefore expect a large variation of 
the measured flux density, whereas above the critical fre- 
quency we expect low flux density variations due to weak 
scintillation. Below /c strong scintillation will occur, with 
a branch due to diffr active scintil lation and one due to re- 
fractive scintillation. iLorimer Kramer! (j2005|) use a dif- 
ferent functional dependence for the transition frequency 
which yiel ds an /c 6 GHz w hich agrees with the es- 
timates by ICamilo et al.l (|2007d ) but is also smaller than 
/c ^ 140 GHz as derived from the NE2001 electron density 
model (jCordes h Lazioll2002! ) . Despite these differences, the 
usage of large integration times and observing bandwidths 
will effectively average a number of scintils, resulting in re- 
liable flux density measurements if band width and integra- 
tion time are sufficiently large. Following Lorimer Kram e^ 
1 20051 ), at 15 GHz we estimate a diffractive scintillation 
bandwidth Azy^z = 320 MHz and a diffractive scintillation 
timescale of id — 600 s which is in good agreement with the 
characteristics of features seen in the dynamic spectrum at 
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the same frequency by iRansom et al ] (|200/t ). Our observa- 
tions always average over several diffr active scintles in fre- 
quency and in time, but will be affected by refractive scintil- 
lations at high frequencies. Around 15 GHz, the modulation 
index for refractive scintillations is estimated as rrir = 0.6 
and the timescale for refractive modulations turns out to 
be tr — 12.8 hours. Our observatories could not track the 
source for such a long time, so that we account for these 
possible variations with an increased quoted flux density 
error. For typical integration times of tint — 40min, our 
observing set-up should result in typical errors of our in- 
dividual flux density measurements due to refractive and 
diffractive scintillations of about 20% at 1.4 GHz, rising to 
a maximum of 77% at 14.6 GHz and levelling off at about 
62% at 33 GHz and above. Similar estimates for the low-DM 
reference source PSR B 1929+ 10 are in good agreement with 
the observations. We therefore conclude that the ISS plays a 
significant role in the individual flux density measurements, 
although it cannot be responsible for the pulse-to-pulse mod- 
ulation index, and the described variations of the spectral 
index for different features of the profile between individual 
sessions. 

Considering the overall flux density of AXP J1810-197 , 
we can divide our observations, spanning more than a year, 
in four intervals. The first lasts until July 2006 when the 
average flux density from the source was above 6 mJy. The 
second lasted from July to September 2006 when the average 
flux density was above 1 mJy. The third epoch lasted from 
October 2006 to July 2007 when the flux density was below 
0.5 mJy and, finally the time after July 2007 when the source 
became too weak for regular detection. At the same time, 
we observe a trend that the flux density variations are larger 
in the low-flux stages. This naturally makes the flux density 
spectrum computation less certain. 

Due to the variability of the flux densities, it is diffi- 
cult to compute an average spectrum. Especially when we 
inspect the flux density measurements from the third epoch, 
we see many cases in which the high and the low frequencies 
vary in a completely different manner for the same day. As 
we have seen, some anti-correlation in flux density variations 
between high and low frequencies appears present, and given 
our discussion above, we consider these variations to be in- 
trinsic to the source and not due to ISS. They can in part 
be explained by the assumption that we observe different 
components of the MP on different days. However, a similar 
variation is also observed for the IP, where always the same 
emission component is visible. Therefore, we conclude that 
in contrast to pulsars, the radio emission of magnetars is 
not intrinsically stable. That is consistent with the overall 
decay of the flux density in recent months and may indicate 
that the radio emission is a transient phenomenon that was 
triggered by the high-energy outburst. As the conditions for 
radio emission may revert back to the pre-outburst stage, 
it will be interesting to monitor the source during and after 
the next outburst. 



5 SUMMARY 

As a result of a coordinated measurement campaign of three 
telescopes operating simultaneously at four different fre- 
quencies and several daily Effelsberg measurements up to 



frequencies of 43 GHz, we flnd AXP J1810— 197 to be an un- 
usual pulsating radio source. Its spectral properties and tem- 
poral fluctuations differ remarkably from normal pulsars. 

A complex picture of the variability of the radio flux 
density emerges as a result of our observations. Signiflcant 
variability exists on all considered time scales, from pulse to 
pulse, day-to-day and over the time of weeks and months, 
most of being it due to intrinsic variations (proflle changes) 
and only some of it affected by scintillations. 

Normal pulsars have stable average proflles enabling us 
to model their emission mechanism. This is however not pos- 
sible with this source. Also the very flat spectrum and the 
visibility to high frequencies makes the source unique among 
radio emitting neutron stars. 
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